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Abstract 

Results from a large sample of hydrodynamical/N-body simulations of galaxy clus- 
ters in a ACDM cosmology are used to simulate cluster X-ray observations. The 
physical modeling of the gas includes radiative cooling, star formation, energy feed- 
back and metal enrichment that follow from supernova explosions. Mock cluster 
samples are constructed grouping simulation data according to a number of con- 
straints which would be satisfied by a data set of X-ray measurements of cluster 
temperatures as expected from Chandra observations. The X-ray spectra from sim- 
ulated clusters are fitted into different energy bands using the XSPEC mekal model. 
The biasing of spectral temperatures with respect to mass- weighted temperatures is 
found to be influenced by two independent processes. The first scale dependency is 
absent in adiabatic runs and is due to cooling, whose efficiency to transform cold 
gas into stars is higher for cool clusters and this in turn implies a strong dependency 
of the spectral versus mass-weighted temperature relation on the cluster mass. The 
second dependency is due to photon emission because of cool gas which is accreted 
during merging events and biases the spectral fits. These events have been quanti- 
fied according to the power ratio method and a robust correlation is found to exist 
between the spectral bias and the amount of cluster substructure. 

The shape of the simulated temperature profiles is not universal and it is steeper 
at the cluster center for cool clusters than for the massive ones. This follows owing 
to the scale dependency introduced by cooling which implies for cool clusters higher 
central temperatures, in scaled units, than for massive clusters. The profiles are in 
good agreement with data in the radial range between ~ 0.1r„j r and ~ 0Ar v i r ; at 
small radii (r ^ 0.1r„j r ) the cooling runs fail to reproduce the shape of the observed 
profiles. The fit is improved if one considers a hierarchical merging scenario in which 
cluster cores can accrete cooler gas through merging with cluster subclumps, though 
the shape of the temperature profiles is modified in a significant way only in the 
regime where the mass of the substructure is a large fraction of the cluster mass. 
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1 Introduction 



Observations show that most of the baryons present in a galaxy cluster are 
in the form of a hot ionized gas (intracluster medium , ICM), with temper- 
atures in the range ~ 10 6 — 10 7 °K (Sarazin 1986). The continuum emission 
is dominated by bremsstrahlung processes, and is free from the contamina- 
tion effects which may arise in the optical band. For this reason, observations 
of galaxy clusters in the X-ray band have been extensively carried out over 
the past decade (Henry et al. 1992, Ebeling et al. 1997, Rosati et al. 1998) 
X-ray observations of the spatial distribution of cluster density and temper- 
ature allow one to determine, under the assumption of hydrostatic equilib- 
rium, the mass in baryons as well as the dark matter profile and the to- 
tal cluster mass Mx- The latter quantity is connected to the global cluster 
temperature Tx, and through the Mx — Tx relation the cluster mass func- 
tion can then be used to find the cluster X-ray temperature or luminosity 
functions. Observations of these quantities are thus used to constrain theo- 
retical models (Henry & Arnoud 1991, Edge et al. 1990, Henry 1997, Rosati 
et al. 1998, Ebeling et al. 1998). Knowledge of the temperature profile is 
also important in order to assess the role of non-gravitational heating pro- 
cesses which can contribute to the budget of the gas thermal energy. This is 
relevant because observations show that the cluster scaling relations do not 
obey the self-similar behavior predicted by gravitational collapse (David et 
al. 1993, Allen & Fabian 1998, Markevitch 1998). Heating of the ICM by 
non-gravitational processes will break self-similarity and help to explain the 
observed relations (Evrard & Henry 1991, Kaiser 1991, White 1991, Loewen- 
stein & Mushotzky 1996, Valageas & Silk 1999, Loewenstein 2001). 

A powerful tool for studying the dependence of cluster observational quantities 
upon theoretical models is given by the use of hydro/N-body simulations. 
Hydrodynamic simulations have been widely used to investigate the formation 
and evolution of galaxy clusters (Katz & White 1993, Suginohara & Ostriker 
1998, Anninos & Norman 1996, Yoshikawa, Jing & Suto 2000, Pearce et al. 
2000, Lewis et al. 2000, Bialek, Evrard & Mohr 2001, Muanwong et al. 2002, 
Loken et al. 2002, Valdarnini 2003, Tornatore et al. 2003, Ascasibar et al. 
2003), see also Voit (2005) and references cited therein. The advantage over 
analytical methods is that they can treat the gas evolution self-consistently. 
The validity of the numerical approach is supported by X-ray observations of 
surface brightness maps, which indicate the existence of merging substructure 
(Mohr et al. 1995). Moreover, high resolution observations performed with 
Chandra and XMM-Newton satellites have revealed the existence of complex 
cluster temperature structure, like shocks induced by mergers and 'cold front' 
(Markevitch et al. 2000, Vikhlinin et al. 2001, Markevitch et al. 2002). 
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Hydrodynamic simulations in which modeling of the gas physical processes 
incorporates radiative cooling, star formation and energy feedback from su- 
pernovae yield a cluster luminosity-temperature relation in good agreement 
with data for cluster temperatures Tx ~ IkeV (Valdarnini 2003, Borgani et 
al. 2004). These results support the radiative cooling model (Bryan 2000), 
where cluster X-ray properties are determined by the cooling efficiency to turn 
cold gas into stars. However, there are a number of issues for which hydro- 
dynamic simulations that include the physical modeling of the gas described 
above do not produce a satisfactory agreement with data. For instance, the ra- 
dial temperature profiles are not isothermal and decline with radius. At small 
distances from the cluster center the simulated profiles have a steep rise mov- 
ing inwards. This behavior is not seen in a set of observations (Markevitch et 
al. 1998, Allen, Schimdt & Fabian 2001a, De Grandi & Molendi 2002, Vikhlinin 
et al. 2005, Piffaretti et al. 2005). The data agree with the simulation results 
at r ^ 20%r vir but indicate a temperature profile which is flat or declining 
at very small distances. These data are also at variance with the predictions 
of the standard cooling flow model (Fabian 1994). In the very central cluster 
regions the amount of cool gas is smaller than that predicted by the model, 
with gas temperatures never below ~ IkeV (Kaastra et al. 2001, Tamura et 
al. 2001, Peterson et al. 2001). This failure to reproduce temperature data has 
prompted many authors to consider possible heating mechanisms to explain 
the lack of cool gas in the cluster core. Among the proposed models, thermal 
conduction (Voigt et al. 2002, Jubelgas, Springel & Dolag 2004, Dolag et al. 
2004, Sijacki & Springel 2006) and energy feedback from active galactic nuclei 
(Churazov et al. 2001, Briiggen & Kaiser 2001, Bohringer et al. 2002, Fabian 
et al. 2002, Briiggen & Kaiser 2002, Omma et al. 2004, Reynolds et al. 2005) 
are the most considered. 

The observed temperature profiles are then sensitive to different physical phe- 
nomena and can be used to constrain theoretical models. Moreover, as already 
outlined, the intracluster gas exhibits a complex thermal structure. It is then 
important to perform properly the comparison of X-ray observations with the 
results of numerical simulations, taking into account the effects of instrumental 
response and background contaminations on the spectral fit temperatures. As 
an additional effect which can bias comparisons, the spectral temperatures are 
measured in a specified energy bands, whereas the simulated cluster tempera- 
tures are theoretically defined according to a specified weighting scheme. This 
problem has been already analyzed by Mathiesen & Evrard (2001) and, more 
recently, by Gardini et al. (2004), Mazzotta et al. (2004), Rasia et al. 2005 
and Vikhlinin (2006). The aim of this paper is to investigate for systematic 
effects which can bias spectroscopic measurements of cluster X-ray tempera- 
tures. For this purpose, a large numerical sample of simulated clusters is used 
to construct spatially resolved X-ray spectra as expected from observations 
with the Chandra satellite. The cluster spectral fit temperatures are found 
using the XSPEC library, and are compared against mean gas temperatures 
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with theoretically defined averages. Moreover, spatially resolved X-ray spec- 
tra of the simulated clusters are also used to investigate how the measured 
temperature profiles differ from the projected profiles obtained directly from 
simulations. The comparisons are performed on a statistical basis, using the 
cluster numerical sample. The numerical sample is also subdivided according 
to the amount of substructure present in a given cluster, this allows us to 
investigate how spectral measurements are affected by the cluster dynamical 
state. 

This paper can be considered a generalization of that of Mathiesen & Evrard 
(2001, hereafter ME), the main differences being the size of the numerical clus- 
ter sample used, together with a more complete physical modeling of the gas 
in the simulations and the X-ray spectroscopic analysis of radial temperature 
profiles. The paper is organized as follows. The simulations are presented in 
Section 2, in Section 3 the procedures with which average and spectroscopic 
temperatures are extracted from the simulated clusters are discussed. Section 
4 is dedicated to the method used to measure substructure in the simulated 
cluster. The results are discussed in Section 5 and the conclusions are drawn 
in Section 6. 



2 Simulations 

The simulations were run using a TREESPH code. A detailed description of 
the simulation procedure is given in Valdarnini (2003, hereafter V03) The cos- 
mological model assumes a flat CDM universe, with matter density parameter 
Q m = 0.3, tt A = 0.7, tt b = 0.019/i- 2 and h = 0.7 is the value of the Hubble 
constant in units of lOOKmsec -1 Mpc~ l . The primeval spectral index of the 
power spectrum n is set to 1 and the power spectrum has been normalized to 
ag = 0.9 on a 8h~ 1 Mpc scale. Initial conditions for the cluster simulations are 
constructed as follows, following a two-step procedure. A collisionless cosmo- 
logical simulation is first run in a box of comoving size L, using a P3M code 
with N p particles, starting from an initial redshift Zi = 10. 

Clusters of galaxies were located at z = using a friends-of-friends (FoF) 
algorithm, so as to detect overdensities in excess of ~ 200f2~ 0,6 within a radius 
r V i r . The corresponding mass M V i r contained within this radius is defined as 
M vir = (in /3)fl m p c A c r^ ir , where A c = 187f2~ a55 for a flat cosmology and 
p c is the critical density. A cluster numerical sample is constructed with the 
cluster identified using this procedure, with the clusters being sorted according 
to the values of their M V i r at z — 0. The final numerical sample used in the 
hydrodynamic simulations consists of two distinct samples, SI and S2. The 
first sample comprises the most massive ri\ = 120 clusters identified at z = 
in a cosmological run with Li = 200h~ l Mpc and N\ p = 84 3 particles. Sample 
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S2 consists of the n 2 most massive clusters identified in a new cosmological 
run with box size L 2 = 400/i _1 Mpc and N 2p = 168 3 . The number n 2 (= 32) of 
clusters in sample S2 is chosen so that the mass M vir of the n 2 —th least massive 
cluster is just above that of the first cluster of the sample SI. For sample S2 
M V i r ranges from ~ 10 15 /i _1 Mq down to 2.3 ■ 10 14 /i~ 1 M Q , while for sample 
SI the nl — th cluster has M v i r = 1O 13 /i _1 M . The random realization of the 
initial density perturbations are different in the two cosmological simulations. 
The final numerical cluster sample of 153 clusters represents then a complete 
sample of the clusters present in a cubic region of size 400/i~ 1 Mpc down to 
virial masses ~ 2.3 • 10 14 /t _1 M Q , and is undersampled by a factor 8 below this 
mass and M vir ~ 10 13 /i _1 M©. 

Hydrodynamic simulations are performed in physical coordinates for each of 
the clusters present in the sample. The initial conditions are constructed as 
follows. All of the cluster particles at z = which are within a sphere of radius 
r V i r located at the cluster center are identified. These particles are traced back 
to a redshift Zi n = 49 and a cube of size L c ~ 25 — 50Mpc enclosing these 
particles is placed at the cluster center. A lattice of N L = 51 3 grid points 
is set inside the cube, and to each node is associated a gas particle with its 
mass and position. A similar lattice is set for dark matter particles, with the 
node positions being displaced by one-half of the grid spacing with respect the 
nodes of the gas lattice. The particle positions are then perturbed, using the 
same random realization as for the cosmological simulations. High-frequency 
modes are added to the original random realization in order to sample the 
increased Nyquist frequency. The gas and dark matter particles used for the 
hydrodynamic simulations are those for which the perturbed particle positions 
lie within a sphere of radius L c /2 from the cluster center. To model external 
gravitational fields these particles are surrounded out to a radius L c by a 
spherical shell of low-resolution dark matter particles, each having a mass 8 
times the sum of the masses of a gas and dark matter particle of the inner 
region. 

The simulations are evolved to the present using a multistep TREESPH code 
with a tolerance parameter 9 = 1, quadrupole corrections enabled, and min- 
imum timesteps of At min = 6.9 • 10 5 yr and At = 8At min for gas and dark 
matter particles, respectively. For the clusters of sample S2(S1), simulations 
have been performed setting e g = 25(15)fcpc, where e g is the gravitational soft- 
ening parameter of the gas particles. The gravitational softening parameters 
of other species of particles have been set according to the scaling e 4 oc m^ 3 , 
where m ; is the mass of the particle i. The runs have a number of gas par- 
ticles N g ~ 70, 000, with similar values for the number of dark matter parti- 
cles in the inner and outer shell. The mass of the gas particles ranges from 
m g ~ 5.5 • 1O 9 M for the first cluster of sample S2, down to m g ~ 6 • 10 8 M Q 
for the least massive cluster of sample SI. For the runs considered here this 
mass resolution can be considered adequate (Valdarnini 2002), yielding final 
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gas distributions with converging profiles for the simulated clusters. This issue 
will be discussed in more detail in sect. 6. 

Numerical integrations have been performed with comoving softenings out to 
z = 20, after which they are kept fixed in physical coordinates. Hydrody- 
namic simulations are followed according to the SPH method (Hernquist & 
Katz 1989), with a minimum smoothing length of s g /A. During its evolution 
the thermal energy equation is subject to energy sinks from recombination 
and collisional excitation, bremsstrahlung and inverse Compton cooling. The 
cooling function A C (T, Z) depends on the gas temperature and metallicity. 
Lock-up tables have been constructed from Sutherland & Dopita (1993) and 
stored in a file. During the simulations, the tabulated values are interpolated 
to obtain A C (T, Z) for the gas particles. Cold gas in high density regions is 
thermally unstable and is subject to star formation (SF). A gas particle which 
is in a collapsing region and whose cooling time is smaller than the dynamical 
time Td is eligible to form in a timestep At a star particle with half of its 
mass, at a rate given by the tJ 1 . At each timestep a star particle heats its 
gas neighbors by supernova (SN) explosions of type II and la. The number of 
SN explosions in the time interval is calculated according to the star particle 
age and the initial mass function, with each SN explosion adding ~ 10 51 ergs 
to the gas thermal energy. An Arimoto & Yoshii (1987) initial mass function 
has been assumed (V03), in the mass range from 0.1 to 40M Q . The gas is also 
metal enriched at each timestep through SN explosions. The mass in metals 
ejected by a star particle is distributed among its neighbors according to the 
SPH smoothing procedure. The gas metallicity thus increases with time, and 
the dependence of the cooling rate on the gas metallicity is taken into ac- 
count properly. In order to construct simulated X-ray spectra from a given 
cluster all the relevant hydrodynamical variables of an individual run are out- 
put during the simulation at various redshifts, with the smaller ones being 
z = 0.025,0.039,0.052. 



3 Temperature definitions 

In this section the way in which mean cluster gas temperatures and tempera- 
ture profiles with theoretically defined averages are obtained from simulations 
data, together with the spectral temperatures, is described. The simulation 
temperatures are defined as those obtained from theoretically defined aver- 
ages applied to the gas particle temperatures of output data. A spectroscopic 
X-ray measurement of a galaxy cluster temperature is constructed with a much 
less simple procedure. Using output data produced with the runs described in 
Section 2, X-ray spectra are first obtained from a certain cluster at a specified 
redshift. As a second step a background component has been added to the 
spectra, and the resulting spectra have been then properly convolved with a 
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template area response file (ARF) and redistribution matrix file (RMF) of the 
Chandra observatory. Finally, the generated event files are fitted by a single 
temperature mekal model available with the XSPEC library. 



3.1 Simulation temperatures 



In the continuum limit an average gas temperature is defined as 
T 1 



nR _ J R T(x)Wd 3 x 

W — r R 



/„" Wd 3 x 



where W is a weight function and the integral has a spherical boundary of 
radius R. Common choices for W are the gas density p g (x) (mass- weighted 
temperatures, W = p g ) and the X-ray emissivity e^ x \x) = A^ X \T, Z)n e nn 
(emission- weighted temperatures, W = £^ x '), n e and uh being the elec- 
tron and hydrogen gas number density, respectively. The cooling function 
AW(T,Z) is calculated using a Raymond-Smith (1977) code. Because of the 
Lagrangian nature of SPH simulations, volume integrals are replaced by a 
summation over the particles, so that 

T£ = EW- r */EW-» (2) 



where the subscript i denotes the value of the quantity at the position of 
the particle i, p = p g , w = m or w = ew according to the chosen weighting 
scheme and the summation is understood to be over all the gas particles 
within a distance R from the cluster center. The cluster center is defined to 
be the position x c where the gas density reaches its maximum. The position 
x c is found according to the following iterative procedure: at each iteration 
k = 1, 2, .., M the center of mass is found, x[ k \ based on the gas particles which 
are within a sphere of radius R^ < i?^ -1 -*, located at x^' 1 ^. The last iteration 
M is such that there must be at least a number Nm of gas particles within 
R( M \ Robust results are obtained for Nm ~ 50 - 100 and i? (fc) = (3/ 'A)R ( - k ~ l \ 
The cluster center is then defined as x\ 

The emission-weighted temperature T e ^ is defined using the bolometric emis- 
sivity, in a similar way a band-limited emission- weighted temperature T^Ei — 
E 2 ) can be defined in the energy range Ei — E 2 by using the band-limited emis- 
sivity e = e u dv, where e v is the specific emissivity. It is useful to define a 
radius i?A such that the average density within that radius is A times the crit- 
ical density, i.e. M A = 4irR 3 A Ap c (z)/3, where p c {z) = 3H(z) 2 /8nG, H(z) 2 = 
HqE(z) 2 and E(z) 2 = Q m (l + z) 3 + Q\. Hereafter, average global cluster tem- 
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peratures will be denoted by 7~, where A takes the values commonly used in 
literature, i.e. A = 2500, 500 and 200. 

Temperature profiles are obtained by first choosing a line of sight and then 
locating in the plane orthogonal to the line of sight annuli of increasing radii 
around the cluster center. Hydrodynamic variables are estimated at a set of 
grid points using the SPH smoothing procedure. The geometry of the grid is 
cylindrical with coordinates (p, <j), z) and the origin is located at the cluster 
center, the z— axis is along the line of sight. The points are uniformly spaced, 
linearly in the angular coordinate, and logarithmically in the radial and z co- 
ordinates. The range of the spatial coordinates is between 2 • 10 -4 and unity 
in units of i?2oo- There are Np = 60, Ni = 20, 2iV 5 = 50 points in the coordi- 
nate intervals. For a continuous distribution the annulus having a finite width 
with projected radial boundaries b\ and 62 has a projected emission- weighted 
temperature T(b) given by 




where b = (bi + &2V2 and \x\ 2 = p 2 + z 2 . The integrals (3) are approximated 
by first evaluating the corresponding hydro variables at the grid points and 
then performing the discrete volume integrals over the set of points which 
are within the given boundaries, the integral along the z— axis is truncated at 
\z\ = -R200 The projected emission- weighted temperatures are calculated at 
a set of 20 projected radii, with approximately uniform intervals in log b and 
ranging between b min ~ 5 ■ 10~ 3 and b max ~ 0.5 in units of i^oo- 

3.2 Spectral temperatures 

The simulated spectroscopic temperatures are obtained by fitting the cluster 
spectral emission with single temperature models. The procedure to obtain 
spectroscopic temperatures from simulated source spectra consists of three 
separate steps which are described in what follows. The first step is to obtain 
photon spectra from a simulated cluster observed at a redshift z along a given 
line of sight. The photon flux per unit energy S u (photons/sec cm 2 keV) that 
reaches the observer is given by 

S u = (l + z) 2 L u{1+z) /A7id 2 L (z), (4) 

where di is the luminosity distance and L v is the volume integral of the specific 
X-ray emissivity of the cluster: L u = J e^dV. The cluster global spectral fit 
temperatures T S A , to be compared with the averaged cluster temperature T A , 
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is obtained from a flux with volume integral performed as in Eq. (2), over 
the same spherical boundaries. The cluster spectral temperature profile T s (b) is 
found using for each ring the photon flux S v (b). The fluxes are calculated with 
the volume integrals over e^> in Eq. (4) being performed in an analogous way 
to those over s^ x > in Eq. (3), using the same grid geometry and boundaries. 

The second step consists of manipulating the simulated photon spectra S v ac- 
cording to a certain procedure, in order to reproduce the spectroscopic temper- 
atures that would be obtained by fitting the cluster spectra as measured by the 
Chandra Advanced Camera for Imaging and Spectroscopy (ACIS) instrument. 
In particular, the spectral analysis described here will consider observations 
relative to the S3 chip of the ACIS configuration. The S3 chip has a small field 
of view (8s3 = 8.4 arcmin), but with a high spatial ( 1 pixel=0.5 arcsec) and 
energy resolution (AE ~ lOOeV^). For the considered energy band the simu- 
lated spectra are then binned into energy channels with width AE = 100e^, 
an energy resolution which is adequate to correctly model the Chandra ACIS 
energy resolution (ME). The cluster global spectral fit temperatures T S A are 
calculated in two energy bands : [2 — 10] and [0.5 — 10]A;eV. Spectroscopic tem- 
perature profiles T s (b) are instead considered in the energy band [0.5 — 7]keV. 
Hereafter it is understood that all the considered energy bands are corrected 
for the cluster redshifts. Once the fluxes S u have been discretized the files 
which contain the spectral distribution S^e are transformed into new files 
Sae{XS), written in the FITS file format 2 . This format allows the files to 
be read by the XSPEC (v 11.3.0) library 3 used to construct the modeled 
spectra. 

The background is modeled by adding a background file B& E (XS) to each 
spectral file Sae(XS). The file Bae(XS) with the proper spectral distri- 
bution of background photons is constructed as follows. The normalization 
parameters are first calculated using, for different energy bands, the quies- 
cent background rates B a (photons/sec arcmin 2 ) taken from Tables 1 and 2 
of Markevitch 4 . For a mock observation with exposure time t exp , aperture 
angle AO and a given energy band, the corresponding background rate is used 
to calculate the expected number of background photons Nb- 

A file Bae(XS) with a background spectrum is then created from a template 
background file 5 using the fakeit command in XSPEC. The spectrum is con- 
volved with the appropriate ARFs and RMFs used to mimic the instrument 
responses (see Sect. 3.4), which are the same pair of matrices used to per- 
form the convolution of the source spectrum Sae(XS). The parameters of 
the background file are rescaled so that when integrating the spectrum over 

2 http://fits.gsfc.nasa.gov/ 

3 htt p : / /heasarc . gsf c . nasa. gov / docs /xanadu /xsp ec/ 
http:/ /cxc. harvard.edu/contrib/maxim/bg/index. html 

5 http:/ /cxc. harvard.edu/cal/Acis/Cal_prods/bkgrnd/current/ 
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t exp seconds there are Nb photons in the energy band under consideration. To 
this end the files are manipulated using the utility fparkey from the FTOOLS 
library 6 . A companion file (XS) is created for later use with the same 
procedure. At this stage the binned source spectrum Sae(XS) is convolved 
with the ARF and RMF files of the mock observation using the atable tool, 
available with XSPEC. The produced pulse height amplitude (PHA) file and 
the corresponding background file are added together by applying to the files 
the mathpha utility of the FTOOLS library. Errors in each energy bin are cal- 
culated using the Poissionian statistic. The resulting spectrum is then grouped 
using the task grppha of the same library, so that the energy channels contain 
at least 20 photon counts. 

Spectral analysis is then finally performed in the third step. The spectrum 
is analyzed with XSPEC, provided that a certain number of constraints are 
satisfied. These constraints take into account the instrumental limits of the 
receiver and are described later (see Sect. 3.3). If the spectrum does not vio- 
late any of these constraints then the fit is performed modeling the emission 
with a single-temperature optically thin plasma. This emission model is the 
single-temperature mekal model implemented with XSPEC and has three free 
parameters: the gas temperature, the metallicity abundance and the normal- 
ization; the other parameters being kept fixed. The spectrum is background 
subtracted using the background file B^ffl previously generated. Errors at the 
68% confidence level (c.l.) for the considered parameters are estimated using 
the tasks error and steppar. 

3. 3 Construction of spectral temperature data sets 

The construction of spectral samples is constrained by the instrumental limits 
of the receiver. For the considered emission models and energy bands the spec- 
tral fit procedure described here is applied to a simulated cluster at redshift z 
if a number of constraints are satisfied. These constraints are chosen so that a 
mock sample of spectral temperatures is generated which satisfies a number of 
properties as those of a set of spectral temperatures obtained from real cluster 
data. These constraints are set as follows. The signal-to-noise ratio (S/N) is 
calculated as 



where P is the source photon rate (photons/sec), B = B a ir(A6) 2 is the photon 
background rate, A9 the source angular aperture and t exp is the observation 
time. For global cluster temperatures the spectral samples are constructed 

6 http:/ /heasarc. gsfc.nasa.gov/lheasoft/ftools/ 




(5) 
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using data from the simulation ensemble at redshifts z = 0.087, 0.47, 0.85. For 
a given energy band the exposure time is held fixed for all the clusters at all the 
redshifts, t exp = 220ks for the energy interval [0.5 — lOj/ceV^ and t exp = 440ks 
for the interval [2 — lOjfceV. In order for the spectral fits to be performed, 
a first criterion that must be satisfied is that the number of source photons 
must exceed a a minimum threshold of JV 7 > 250. Another criterion is that the 
cluster spectra must have a signal-to-noise ratio of at least S/N > 10. This 
threshold selects spectra of good quality, and at high redshifts this prevents 
the spectral fitting of background dominated spectra. 



Spectral fits are performed in the energy band [2 — 10] keV only for those 
cluster with T^® > 2keV. Spectral temperature profiles are constructed from 
simulation data at redshifts z = 0.116,0.052,0.039 and z = 0.025. For these 
spectra t exp = 140/cs for all the clusters and redshifts in the considered energy 
band [0.5— 7]keV. The spectrum of each annulus is fitted if it satisfies iV 7 > 250 
and B < 2P. The last criterion extracts from the ring under consideration 
spectra having S/N with high values (~ 10 2 ). This removes from the spectral 
sample fits of poor quality, though a large number of spectra can still be fitted. 
For global cluster temperatures the redshifts have been chosen in order to 
highlight the biases associated with the cluster selection owing to the detector 
physical limits. Therefore there are three different and well separated redshifts 
rather than the redshift distribution which would be expected from a real 
cluster catalog. For spectral temperature profiles the adopted range of redshifts 
was motivated by the need to mimic that of the clusters for which temperature 
profiles are measured (see, for example, Vikhlinin et al. 2005). 



The following constraints take into account the finite field of view of the S3 
chip. For global cluster temperatures spectral analysis is performed only if 
#a = R^/dA^z) < 9 S3 /2 = 4.2', here (Ia{z) is the angular distance in a flat 
cosmology. For the spectral temperature profiles a ring spectrum is fitted only 
if the ring has inner radius Ti n /dA(z) > 4" or outer radius r out /dA(z) < 4', and 
thickness (r out — r in )/dA(z) > 5". Finally, a spectral fit temperature is rejected 
if, for the considered number of degrees of freedom (d.o.f.), the statistic of the 
fitted model a gives a value of the x 2 larger than the threshold value x\ f° r 
which the null hypothesis probability is larger than 10 %. Values of x\ were 
pre-calculated as a function of the number of d.o.f. up to d.o.f. = 1,000 and 
stored in a table which was used at run time. The range of values for the 
reduced xl = X 2 / d.o.f. of the fits lies between 0.5 and unity, with the bulk of 
the values clustered around 0.5. 
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3.4 Construction of the Response Functions 

The instrumental response of the ACIS-S3 detector is position dependent, but 
can be considered approximately constant within each of the 32x32 square 
regions covering the chip. For extended sources, weighted ARFs and RMFs 
are then constructed by summing the corresponding subregion matrices, with 
a weight proportional to the photon count of the spectral image in the re- 
gion under consideration (Gardini et al. 2004). Here a simplified treatment is 
adopted, in which the source emission is first integrated spatially, and then 
the single source spectrum is convolved with appropriate template ARFs and 
RMFs. For an isothermal gas the two methods give the same results, if there 
is a temperature gradient then photons far from the cluster center will be in 
proportion weighted more than in the procedure which sums over subregions. 
Therefore, this implies on average global cluster temperatures biased toward 
lower values, with respect the first method. Because the bulk of the emission is 
concentrated at the cluster center, this effect is however expected to be negli- 
gible. The same arguments apply to the photon spectra extracted from cluster 
annuli, for which it is assumed that spatial variations of the gas temperature 
within a single annulus can be neglected. 

If one considers a large number of spectra, the weighting errors associated 
with a single convolution can be statistically reduced if the source spectra 
are convolved with weighted ARFs and RMFs extracted from a template ob- 
servation of a cluster with a regular gas distribution. This cluster must have 
the characteristic of being in a highly relaxed state and of covering a large 
field of view in the ACIS-S3 chip. To this end Abell 2029 was chosen, for 
which analysis of the ACIS-S3 image (Lewis, Stocke & Buote 2002) shows 
a very regular structure and a spatial extent ~ 4' wide. The cluster is lo- 
cated at z — 0.0767, spectral analysis yields an intracluster gas temperature 
Tx — 8keV and a virial radius r vir ~ 2.5Mpc. The cluster exhibits a regular 
X-ray morphology and there is no evidence of a cooling flow at the cluster 
center. A set of weighted ARFs and RMFs was generated with different geo- 
metrical boundaries using the tool acisspec of the CIAO 7 version 3.0.2, with 
calibration database CALDB (version 2.28). The selected source regions are 
centered on the peak of the X-ray emission, with their geometry being chosen 
to approximately fit that of the analyzed spectra. For each source spectrum, 
the pair of response matrices with geometrical boundaries providing the best 
fit to those of the considered source is chosen in order to be convolved with 
the binned source spectrum S&e(XS). The same pair is later used when the 
spectrum is fitted using the single-temperature emission model mekal imple- 
mented with XSPEC. The main advantage of the procedure described here 
is that the adopted approximation allows one to analyze with a reasonable 
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computational cost a large number of spectra, as it is in this paper. 



4 Statistic of cluster substructure 

Spectroscopic measurements of the cluster X-ray temperatures can also be 
affected by the cluster dynamical state (ME). The amount of substructure 
present in the inner mass distribution of a galaxy cluster is closely related 
to its dynamical state (Richstone, Loeb & Turner 1992) and many statistical 
measures have been proposed to quantify cluster substructure (Buote 2002 , 
and references cited therein). Analysis of X-ray images does not suffer from 
projection effects which are present in the optical band and the X-ray surface 
brightness is expected to qualitatively follow the morphology of the projected 
cluster mass density (Buote 2002). For this reason, in this work is adopted 
the power ratio method (Buote & Tsai 1995) as a statistical indicator of the 
cluster dynamical state. The method has been widely used to study cluster 
X-ray morphologies (cf. Buote 2000). An application of this method to analyze 
global morphologies of simulated clusters was already performed in an earlier 
paper (Valdarnini, Ghizzardi & Bonometto 1999), and a systematic statistical 
analysis of the evolution of cluster X-ray morphology using hydrodynamic 
simulations in different cosmological models is in preparation (Buote et al. 
2006). 

The method works as follows. The X-ray surface brightness Ex(p, <p) along a 
given line of sight is the source term of the pseudo potential ^(p,ip) which 
satisfies the 2-D Poisson equation. The pseudo potential is expanded into plane 
harmonics and the m — th coefficients of the expansion are given by : 

a m = J d 2 x Ex(f )R cos(m<£> ), (6) 
r'<r 

Pm= J <i 2 x'S x (f')i?' m sin(m<p'), (7) 
r'<r 

where x = (p, ip) and the integration is over a circular aperture of radius R, 
which is also termed the aperture radius R = R ap . The m — th power ratio is 
then defined as 

rt m \R ap ) = \og 10 (P m /P ) , (8) 
where 
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P m (Ra P ) = ^(a 2 m + P 2 J m>0, (9) 
P = [a ln(R ap /kpc)] 2 . (10) 

The ratio P m /Po is a measure of the amount of structure present on the scale 
of the aperture radius R ap . For a relaxed configuration lT m ) — > — oo. The 
values of P m depend on the choice of the coordinate system. If the origin is 
the center of mass then P\ vanishes and is a measure of the degree of 
flattening (Buote 2002). Large values of P3/P0 indicate asymmetric distribu- 
tions and this ratio will be used in the subsequent analysis as an indicator of 
the amount of substructure present in a cluster. The moments P m {Rap) of a 
chosen cluster at a given redshift are calculated along a line of sight as in Eq. 
2, by performing the integrals (6), (7) according to the SPH prescription. The 
origin of coordinates is set at the peak of the X-ray emission. 

In order to analyze how spectral fits are affected by the amount of substruc- 
ture present in a cluster, three separate cluster sub-samples were extracted 
from the cluster sample. The sub-samples were generated by first construct- 
ing the cumulative distribution of the sample values of II3, and then identify- 
ing those clusters with a value of II3 below the threshold values which define 
the 25%, 50% and 75% percentile of the cumulative distribution, respectively. 
Cluster sub-samples will be cross-correlated with values of the spectroscopic 
temperatures in order to investigate how the cluster dynamical state can bias 
spectroscopic measurements. 



5 Results and discussion 

This section is dedicated to analyzing the dependence of spectral fit temper- 
atures, as obtained from simulation data using the procedures described in 
sect. 3, against cluster temperatures defined according to weighted averages. 
Global cluster temperatures are discussed in the first two parts, projected 
temperature profiles from individual clusters are presented in the third part. 

5.1 Global cluster temperatures: effects of cooling 

The spectral fits temperatures T S A are obtained by applying the prescriptions 
of sect. 3 to a cluster sample which is constructed by grouping data from the 
simulation ensemble at redshifts z = 0.087, 0.47 and z = 0.85. Those clusters 
for which the values of II3 exceed the threshold value which defines the 50% 
percentile of the cumulative distribution have been removed from the cluster 
sample. For those clusters which are part of a cluster sub-sample, identified 
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by a given value of redshift z and overdensity A, the cluster power ratios were 
calculated according to the procedures described in sect. 4, by choosing the 
values of the aperture radius R ap to match the radius ( see sect. 5.2). The 
power ratio cumulative distributions were then constructed for these cluster 
sub-samples. This procedure has the effect of removing from the cluster sample 
most of those clusters undergoing a major merger event, thus disentangling 
the effects of cooling from those due to merging when analyzing the biasing 
of spectral fit temperatures . The choice of the threshold value is somewhat 
arbitrary, nonetheless it has been found that the results obtained are fairly 
robust against changes in the chosen value. 

As far as concerns the dependence of the sample on redshift, at high redshifts 
the construction of the sample is mainly limited by noise, for example the 
most massive cluster at z = 0.85 has r 2 oo — l-2Mpc and in the [0.5 — 10]keV 
energy band it emits iV 7 ~ 4, 200 photons. The photon count rate of the source 
is c/s ~ 2 ■ 10 _2 sec _1 , with a similar value for the background rate. The latter 
depends on the source size through B oc A6 2 . At low redshifts the sample 
construction is limited by the geometrical limits and the main contribution 
is given by the low temperature clusters. In the [2 — 10]keV bandpass the 
sample construction is mainly constrained by the photon threshold JV 7 > 250, 
because of the low number of photons emitted by the clusters. In the following 
analysis mass-weighted temperatures, which are supposed to be fair tracers 
of the virial values, will be taken as reference temperatures with which to 
correlate emission- weighted temperatures or spectral temperatures. For a fixed 
cluster overdensity A linear fits of logarithmic variables will be of the form 
log 10 s = a + 61og 10 T£ w , with errors on spectral fit temperatures given by 
the fit of the XSPEC mekal model and those on the best fit parameters will 
be assumed at the 68% c.L. 

The correlation between emission-weighted temperatures and spectral temper- 
atures versus mass-weighted temperatures are displayed in Fig. 1 for different 
overdensities and energy bands. Within a different panel different symbols are 
for different overdensities A, the lines associated with the A symbols corre- 
spond to the linear fits, with the coefficients a and b given in Table 1. For the 
sake of clarity random subsamples are plotted. There is a robust correlation 
between and T^ w with T 2 ^ 00 ~ T^®°. At smaller overdensities this is not 
satisfied, with T 2 °° being systematically higher than T™. Emission- weighted 
temperatures weight preferentially the inner core regions and change a little 
when going from A = 2500 to A = 200, whereas drops to smaller values 
when R A gets higher. From Table 1 T 2 °° /T 2 ™ ~ 1.35(T 2 ™)-°- 01 , so that T 2 °° is 
~ 30% higher than T^, with a weak scale dependence on T™. These results 
are at variance with those of ME. A comparison with their Fig. 5 shows that 
the relation T ew versus T mw is approximately around the line of equality, with 
some preference for the scattered points to lie above the line, i.e. T mw £ T ew . 
These differences must be ascribed to the different physical modeling of the 
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Fig. 1. Spectral temperatures T s are plotted as a function of the cluster 
mass- weighted (T mw ) or emission- weighted (T ew ) temperatures. Upper panel (a) 
refers to T ew versus T mw . For T s the numbers in parenthesis indicate the energy 
band in keV. The temperatures are found according to the procedures described 
in Sect. 3 and are defined within cluster overdensities A = 2500 (filled triangles), 
A = 500 (open circles) and A = 200 (filled squares). The diagonal line is the line of 
equality, the line associated with the A symbol corresponds to the linear fit given 
in Table 1. For the sake of clarity not all of the points and not all of the error bars 
are shown. 

gas in the hydro simulations. In ME the gas was treated adiabatically, here the 
runs incorporate radiative cooling and star formation, as well as SN feedback. 
In the cluster central region the gas is subject to radiative losses and cold 
gas which is Jeans unstable is removed because of star formation. Previous 
analyses (Lewis et al. 2000, Valdarnini 2002) have shown that taking into ac- 
count these effects has the consequences of changing the overall temperature 
profile, with a temperature decline with the radius being steeper than in the 
adiabatic runs. This is the main reason for the discrepancy found between 
the behavior of T mw versus A here and in ME. It must be stressed that this 
discussion refers to relative differences between mass and emission-weighted 
temperatures, when compared against the corresponding ones in the adiabatic 
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runs. 



It is useful to perform a comparison between absolute values. The relation 
between mass-weighted temperatures of different runs can be calibrated using 
the mass-temperature relations. For A = 500 Eq. (6) of ME gives M 500 — T^®. 
The analogous relation here is ( for h=l) 

log 10 E(z)M 500 = (13.39 ± 3 ■ 10" 4 ) + (1.68 ± 8 ■ 10" 4 ) \og w T™. (11) 



To ease the notation in the following discussion the superscript A is dropped 
from the temperatures, the letter c stands for the cooling runs performed 
here and a refers to the adiabatic runs of ME. By equating the two mass- 
temperature relations one finds 

T mw {a) = (0.77 ± 0.04)T miu (c) L06±a02 , (12) 



this allows us to relate T ew (a) to T mw (c) using the linear fits coefficients of 
Fig. 6 of ME. This yields 



T ew (a)/T ew (c) * (0.52 ± Om)T mw (c) om±om . (13) 



So that cooling runs have emission-weighted temperatures which are for A = 
500 a factor ~ 2 higher than in the corresponding adiabatic runs. Similarly, 
the ratio for mass-weighted temperatures is 

T mw {a)/T mw {c) ~ (0.77 ± 0.04)T mw (c) - 07±a01 . (14) 



Therefore, mass- weighted temperatures are then higher by a factor ~ 1.4. 
The ratio T ew /T mw is higher than in ME by a factor ~ 30%, with a weak 
dependence on T mw . It is found 



T 



(a) 



I 



T 



0.67T„ 



\0.125±0.01 



(c) 



(15) 



These relations show how incorporating into the simulations a more realistic 
physical modeling of the gas leads to differences in the gas temperature and 
distribution which can have a strong impact on temperature averages. This 
suggests that also the relationships between spectral temperatures and mean 
gas temperatures can be modified when the simulations take into account 
the effects of radiative cooling and star formation. The relationships between 
spectral fit temperatures and averaged gas temperatures are shown in the 
other panels of Fig. 1. 
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Table 1 

Values of the best-ht coefficients and la confidence limits for linear fits of the form 
logT^ = a + blogT^ applied to the cluster emission- weighted, mass- weighted and 
spectral temperatures in keV units. The coefficients are calculated for cluster tem- 
peratures defined within radii enclosing cluster overdensities A = 2500, 500, 200. 
The fits are performed over samples which are constructed grouping data from the 
simulation ensemble at redshifts z = 0.087, 0.47, 0.85. The subscript w = ew, mw, s 
indicates the weighting schemes, described in Sect. 3. For spectral temperatures the 
numbers in square brackets indicate the energy band. An asterisk in the column of 
the linear fit coefficient a means that the x 2 probability of the fit is below 0.01. 
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relationship between T s and T mw for A = 2500 has many points which are 
approximately distributed according to the linear regression fit, nonetheless 
the x 2 probability of the fit is below 0.01. This follows because, as a conse- 
quence of selection effects, there are many points of good statistical quality 
which are at low temperatures. The best-fit parameters are mostly weighted 
by these points, and the best-fit line follows their distribution. This is more 
clearly illustrated in Fig. 2, where in the left panel (a) are shown all of the 
points with which the spectral sample T s [0.5 — 10] — T mw for A = 2500 is 
constructed. Only spectral temperatures with x 2 values from the spectral fits 
below the threshold value x 2 have been plotted. For a given number of d.o.f. 
xl is such that the x 2 probability of the fit gives P x 2 (d.o.f.) = 10%. As pre- 
viously outlined at high z's the construction of the sample is constrained by 
noise, with only massive clusters being selected, whereas at low redshifts cool 
clusters are preferentially selected. The poor quality of the linear fit is mainly 
determined by the relatively large number of outliers which are present at high 
z's in the spectral sample and which do not follow the distribution of the low 
temperature clusters. This feature is closely related to the way in which the 
introduction of cooling affects the gas distribution. 

ME argue that spectroscopic temperatures T s [0.5 — 10] are biased toward 
lower values of the mass-weighted temperatures because spectroscopically de- 
termined temperatures are weighted by the fitting process according to the 
photon counts, which are dominated by the low-energy part of the spectrum. 
From the spectral sample of Fig. 2(a) it can be seen that for massive clusters 
spectral temperatures are lower than mass-weighted temperatures, whereas 
T s ~ T mw as cool clusters are considered. This behavior follows from two ef- 
fects: the way in which the sample is constructed and the introduction of the 
physical modeling of cooling in the simulations. At high redshifts the sample is 
dominated by massive clusters for which T s £ T mw , as in ME, at low redshifts 
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Fig. 2. In panel (a) it is shown for A = 2500 the relation between T s [0.5 — 10] and 
Tmw° °f Fig. 1(d), different symbols refer to clusters selected at different redshifts. 
The dashed line is the linear fit of Table 1. (b): The source photon distribution 
emitted within i?2500 at z = 0.473 is binned as a function of the photon energy in 
units of the average gas temperature. Continuous (dashed) line is the average over 
the 20 most (least) massive clusters of the sample. Vertical units are arbitrary. 



the sample population is dominated by cool clusters. For these clusters, accord- 
ing to the cooling scenario (Bryan 2000, Voit & Bryan 2001), the efficiency of 
galaxy formation is higher than in hot clusters. This implies a removal of the 
low-entropy cooled gas, transformed into stars, and a subsequent inflow of the 
surrounding high-entropy gas. Therefore, in this scenario, it follows that for 
cool clusters the central cluster temperature, in units of a characteristic cluster 
temperature, is higher than for massive clusters. This issue will be discussed 
in more detail in the section dedicated to analyzing the temperature profiles. 
Fig 2(b) shows the energy distribution of the photons emitted at z — 0.473 
within -R2500 by the gas particles of the the 20 most (least) massive clusters of 
the spectral sample. Energy is in units of the averaged gas temperature within 
-R2500 : < T2500 >— J2iTi/N, where the summation is over all of the N gas 
particles within the radius. 

The histograms of the two distributions clearly indicate that cool clusters will 
have average temperatures higher, when rescaled to a characteristic mean gas 
cluster temperature, than the corresponding ones of massive clusters. Spec- 
troscopic temperatures too are expected to follow this behavior and, as a 
consequence, will have a dependence on mass-weighted temperatures different 
from that found by ME. The argument can be quantified in more detail by 
applying again the same arguments which lead to Eq. 12. From Fig. 6 of ME 
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it is found in the [0.5 — 10] bandpass for A = 500 



T s (a) = (0.6 ± 0.05)T mu) (c) 



(1.1±0.03) 



(16) 



and 



T s (a)/T s (c) ~ (0.56±0.05)T mw (c)(°- 13±0 - 06 ). 



(17) 



Using this equation T s (c) ~ 1.8T s (a) for T mw = IkeV but T s (c) ~ 1.3T s (a) 
at T mw = lOkeV. This shows that incorporating cooling into the simulations 
has the net effect of introducing for spectral temperatures a scale dependency 
much stronger than in the adiabatic runs of ME. The runs performed here 
take into account the dependence of the cooling function on the gas metal- 
licity; line emission from cold gas is expected to bias significantly spectral 
temperatures toward lower values (ME), from the results found here it turns 
out that this bias is largely covered by the scale dependency introduced by 
cooling. The scale dependency of the T s — T mw relationship is also dependent 
on the energy bandpass, as it is found in ME, when passing from the energy 
interval [0.5 — 10]keV to [2 — 10]keV panel (b) of Fig. 1 shows that cool clusters 
have now T s £ T mw at A = 2500. This follows because the low energy photons 
which weight spectral temperatures towards lower values are removed when 
the spectral fitting is performed in the [2 — 10]keV range. 

These combined effects on the scale dependency of T s versus T mw are equally 
present when the cluster radius -Ra, within which the temperatures are found, 
is increased. Here the scale dependency is also affected by the biases associated 
with the cluster selections because of the geometrical constraints that limit 
the sample construction. Thus, the cluster sample with A = 200 of Fig. 1(d) 
is expected to be less populated by massive clusters than the corresponding 
sample with A = 2500. This is because most of the clusters at z — 0.087 
have #2oo greater than the field of view of the S3 chip. This bias is relatively 
unimportant in the determination of the slope of the T s — T mw relationship 
for A = 200, which is dominated by the general tendency of having spectral 
temperatures higher than mass-weighted temperatures. This behavior was al- 
ready detected when discussing the relationship between T ew and T mw versus 
A and it is a consequence of incorporating cooling into the simulations. Ac- 
cording to the simulation results mass-weighted temperatures have a strong 
dependency on A, whereas T s or T ew remain relatively unaffected when 
is increased because their weighting scheme favors the central regions where 
the bulk of the emission is located. Therefore it turns out that spectral tem- 
peratures are systematically higher than mass-weighted temperatures when 
A ;$ 10 3 . These results are at variance with those of ME, for which spectral 
fit temperatures are biased against mass-weighted temperatures by a ~ 20% 
toward lower values. 
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There is also a well defined tendency to have large statistical errors on the spec- 
tral fit temperatures when the overdensity A decreases. This effect is present 
in all of the panels of Fig. 1, where a spectral T s is correlated with a weighted 
average gas temperature. The errors a a and on the best fit parameters are 
given in Table 1 and indicate this tendency. These large errors on the spectral 
T s as Ra is increased are likely to be induced by the inclusion within the 
spectral fit volume of an increasing amount of gas at temperatures lower than 
those at the cluster center. In the [2 — 10]keV bandpass the cluster spectral 
temperatures with A = 200 (Fig. 1(b)) have large statistical errors. The slope 
of the T s — T mw relationship is greater than one, but with a statistical error 
as large as the slope itself. As a general rule it is found that in the energy 
range [2 — 10]keV, spectral fits temperatures have large statistical errors. This 
is mainly due to the poor quality of the photon statistics associated with the 
fits. This is closely related to the construction of the spectral temperature 
data sets. There is an important difference between the way in which spectral 
temperature data sets are constructed here and in ME. Each spectrum is fitted 
by ME, using an isothermal spectral fit model mekal with fixed metallicity, 
keeping the number of observed photons constant to 20, 000. At variance with 
ME, here the exposure time is kept fixed, so that the observed number of 
photons can vary from cluster to cluster. This choice was motivated by the 
request of creating spectral temperature sets from fake observations with re- 
alistic exposure times, therefore avoiding the selection biases associated with 
the inclusions in the spectra data set of cluster spectra with an otherwise low 
photon count. As a consequence, spectral fits in the [2 — 10]keV energy band 
have generically a poor statistic, because of the smaller number of photons 
with respect to the [0.5 — 10]keV band, even with the very high exposure time 
(see sect. 3.3) of 440ks. 

In Fig. 1(c) for different overdensities A the relationships between T s [0.5 — 
10] and T ew are displayed. These have the same scale dependencies as those 
between T s and T mw . The correlations between T s [0.5 — 10] and T ew can then 
be explained using the arguments previously discussed. As a general result it 
appears that T s -T ew relationships are more widely scattered around the line 
of equality than their counterparts T s -T mw . This effect is particularly severe 
when the overdensity A = 200 is considered. This is most likely due to the 
combined effects of the inclusion within the spectral fit volume of clumps of 
cold dense gas as increases and to the p 2 dependence of the T ew measure. 
There is a tendency for the ratio of T s to T ew to be higher than that of T s 
to T mw , because of the biasing between T ew and T mw . A comparison with 
previous results can be performed for the relation T s [0.5 — 10]-T eiu at A = 500. 
From a cluster sample obtained from a set of hydrodynamical simulations, 
with a physical treatment of the gas similar to that used here, Rasia et al. 
(2005) found a relationship at A = 500 between T ew and a suitably averaged 
temperature T s i. The latter is supposed to be a good approximation to the 
spectroscopic temperature T s [0.5 — 10]. A linear fit of the form T s i = aT ew + b 
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Fig. 3. Relative differences between spectral or emission-weighted temperatures ver- 
sus mass- weighted temperatures. The temperatures are found as in Fig. 1 and differ- 
ent symbols correspond to different cluster overdensities : A = 2500 (filled triangles), 
A = 500 (open circles) and A = 200 (filled squares). The notation 5{y)/x stands 
for (y — x)/x. As in Fig. 1, random subsamples are plotted to avoid overcrowding. 



yields a = 0.7±0.01 and b = 0.29±0.05. This relationship is in good agreement 
with the corresponding one in Table 1. At T ew = IkeV the values of T s [0.5 — 10] 
and T s i are nearly equal, whereas at T ew = lOkeV there is a ~ 10% difference 
between the two values. 



The scale dependencies of the quantities shown in Fig. 1 are more clearly 
illustrated in Fig. 3, where the fractional differences between correlated quan- 
tities are plotted. Panel (a) shows that the quantity 5(Ts(0.5 — 10))/T mw for 
A = 2500 is characterized by negative values for T mw ^ 3keV, and by a plume 
of positive values as long as T mw — > IkeV. This feature is exacerbated when 
considering spectral temperatures in the energy band [2 — lOjfceV" ( panel (b) 
), as expected because of the lack of low-energy photons. These results can 
be compared with Fig. 8 of ME, for which the condition T$ < T mw is always 
valid in the [0.5 — lOjfceV^ bandpass. The behavior of 5(Ts(0.5 — 10)/T mw for 
A = 2500 is reproduced also when considering temperatures defined within 
overdensities A = 500 or A = 200, here the distributions are characterized 
by a smaller number of massive clusters because of selection constraints. It 
is clear that the scale dependency of spectral temperatures against mass- 
weighted temperatures is dominated by the effect of introducing cooling in 
the simulations. The biasing of spectral temperatures associated with line 
emission from cool gas was dominant in previous runs which considered only 
gas shock-heating, but it appears now as a secondary effect. 
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5.2 Global cluster temperatures: effects of substructure 

The discussion of the relationship between spectral and mass averaged tem- 
peratures has not yet considered the role of merger events during the individ- 
ual cluster formation history. According to ME accretion and merger events 
will alter the measurement of spectral temperatures. These are expected to 
be significantly biased during a major merger event. As the gas gets shock- 
heated because of the collision its temperature and luminosity will temporarily 
increase. The merging will proceed with the subclump of cold gas moving to- 
ward the cluster center and with its presence will bias the cluster spectrum 
toward lower temperatures. The gas of the clump is expected to be cooler than 
most of the gas in the cluster because of its smaller mass with respect to that 
of the cluster. ME argue that spectral temperatures will be much smaller than 
mass- weighted temperatures in the proximity of a merger event. 

The substructure statistic used here is the power ratio method described in 
sect. 4. This statistic is different from the one used by ME and has the main 
advantage of giving a measurement of the amount of structure present on a 
given scale, chosen accordingly to the value of the aperture radius R ap . This 
allows us to correlate, for different overdensities or radii the fractional 
differences 5(T mw )/Ts as a function of the power ratio n 3 = log 10 P 3 /Po- For 
the clusters under consideration, the values of R ap have been chosen to match 
that of the radius Pa- The two panels of Fig. 4 refer to S(T mw )/Ts as a func- 
tion of n 3 for Ts in the energy bands [0.5 — 10] and [2 — lOjfceV. In each panel 
the distributions as obtained by considering the three different overdensities 
previously considered are displayed. A striking result of Fig. 4(a) is the robust 
correlation which is found between S(T mw )/Ts and the substructure as mea- 
sured by n 3 . A result which is likely due to the theoretical framework of the 
statistical method used to measure substructure, which has the advantage of 
introducing well defined scale dependent quantities. In the [2 — lOjfceV band- 
pass there is a significant correlation only for A = 2500 and A = 500 , the 
other distribution (A = 200) being scattered similarly to the corresponding 
one of Fig. 3(b). 

The correlations are confirmed from the values of Table 2 in which are reported 
the Spearman rank correlations and significance levels of the corresponding 
panels. The correlation behavior reveals that Ts is significantly lower than 
mass-weighted temperatures as n 3 approaches zero, i.e. whenever the cluster 
gas distribution is strongly perturbed, which is an indicator of a merging event. 
As long as Il 3 becomes very large and negative the cluster is in a relaxed state 
and spectral temperatures Ts are higher than mass-weighted temperatures. 
For A = 200 two clusters have been chosen as representative of these two 
regimes, they are identified by the circle drawn around their points; their 
phase-space diagrams are illustrated in Fig. 5. 
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Fig. 4. Relative differences between mass-weighted temperatures and spectral tem- 
peratures are shown versus the cluster power ratios H3 = logioP?,/ Pq. The temper- 
atures are defined as in Fig. 1 within different cluster overdensities A. The values 
of II3 are calculated at the same redshift and line of sight of the corresponding clus- 
ter temperatures. The aperture radii are defined according to the value of A, with 
R ap = -R20O] Rap = -R200/2 and R ap = i?20o/4 chosen to approximately match the 
value of i?A in correspondence of A = 200, 500, 2500, respectively. The two clusters 
identified by a circle in panel (a) have their phase space diagram plotted in Fig. 5. 

The phase-space diagrams of the two clusters in a different dynamical state 
are displayed in Fig. 5. Gas particle temperatures Tj are in units of and 
the solid line gives Tj = 2keV. Left panel (a) is for the cluster identified 
by the circle lying at the bottom left corner of Fig. 4 (a). The cluster is in 
a fairly relaxed state with = 1.8keV, II3 = —7.75 and will be termed 
as 'quiescent'. In the right panel (b) the logT — logra points for the cluster 
marked by the circle in the top right part of Fig. 4(a) are plotted. For this 
cluster T™ = 5.8keV and IT3 = -4.52. The cluster can be considered strongly 
asymmetric along the chosen line of sight, its value of IT3 being above the 
threshold value denning the 75% percentile of the cumulative distribution. 
This cluster will be denoted as 'active'. The main difference between the two 
distributions is the long tail of cool gas which characterizes the active cluster, 
for which T s ~ 0.6T^. Note that the quiescent cluster has most of the material 
at temperatures below the solid line indicating 2keV, nonetheless for this 
cluster T s ~ 

These findings support the analysis of ME, for whom clusters with an under- 
going merging activity have spectral temperatures smaller than mass-weighted 
temperatures. It must be stressed that the scale dependencies of spectral tem- 
peratures which are shown in the panels of Fig. 4 are independent from those 
function of mass-weighted temperatures which follow by introducing cooling 
in the simulations. According to these results, the biasing of spectral tem- 
peratures can be described in terms of a two-parameter model. A first scale 
dependency is introduced by cooling and is a function of the cluster mass 
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Table 2 

Spearman linear correlation coefficients r s and significance levels P ra are given for 
the fractional temperature differences displayed in the two panels of Fig. 4 as a 
function of the cluster power ratios H3. Values of P rs below 10 -2 have been rounded 
to zero and a significance level of 5% is used to reject the null hypothesis r s = 0. 
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Fig. 5. Phase-space diagrams of the two clusters identified in Fig. 4(a) by the circles 
drawn around the plotted points. Left panel (a) is for the cluster lying at the bottom 
left part of the panel. Temperatures are in units of T^® = 1.8 (5.8)keV for the left 
(right) cluster. The number density ratios are defined as rii/n c = pi/p c (z). Filled 
symbols are for gas particles within -R2500 and open symbols for those within i?200- 
The continuous line gives 7$ = 2keV. 

or mass-weighted temperature. The second dependency is correlated with the 
amount of cluster substructure and is independent from the first. Given the 
importance of this argument it has been decided to look at the phase-space 
diagram of a cluster with its position in Fig. 4(a) as close as possible to the 
one identified by the active cluster, but characterized by a low value of T^. 
A cluster was chosen with T s ~ 0.7T^™(~ IkeV) and n 3 = -4.36. For this 
cluster, a distribution of points in the log T — log n plane very similar to that 
of the active cluster was found. This confirms that the two processes which 
govern the biasing of spectral temperatures, subclump accretion and radiative 
cooling, introduce scale dependencies which can be considered independent of 
each other. 



25 



Finally, there is a further source of biasing which is due to selection constraints 
because of the detector configuration. At high redshifts the sample construc- 
tion is limited by noise, whereas at low redshifts the detector geometry is 
dominant. The relationships between spectral and mass- weighted tempera- 
tures are then shaped by these selection effects, which weight the scale depen- 
dency introduced by cooling according to the cluster redshift. Spectroscopic 
measurements with other X-ray telescopes will give relationships different from 
those obtained here using the Chandra ACIS-S3 detector, though qualitatively 
similar. Because of the limits to the instrument sensitivity the generic popu- 
lation of a cluster sample will be dominated at low redshifts by cool clusters 
and at high redshifts by massive clusters. Another issue concerns the geomet- 
rical limits of the detector, if its field of view is larger than that of the chip 
S3 spectral relationships for A = 200 are expected to include a number of 
massive clusters larger than the one of Fig. 1 (d). This in turn implies a scale 
effect which will change the bias between T$ and T mw . This effect is however 
expected to be negligible since most of the new clusters will be characterized 
by T s ~ T mw . Thus the scale dependencies of spectral temperatures found here 
are expected to be qualitatively reproduced with other detectors. 



5.3 Radial temperature profiles 



In this section the behavior of the radial temperature profiles as obtained from 
the sets of simulated spectral samples constructed according to the procedures 
described in Sect. 3 are discussed. Spectral temperature profiles T$(r) are com- 
pared against projected emission- weighted temperature profiles T ew (r), the 
spectral fits are performed in the energy band [0.5 — 7]keV and the emission- 
weighted temperatures are defined in the same energy band. The emission- 
weighted profiles have been calculated keeping fixed the chosen line of sight 
and the radial binning, in units of r2oo, for all of the clusters. For a given 
cluster, the corresponding spectral temperature profile is obtained by fitting 
the photon spectra of the considered rings. Because of the scale dependencies 
previously discussed, the sample of temperature profiles has been subdivided 
by grouping individual profiles into sub-sample according to several cluster 
properties. A cluster is part of a subsample denoted by 'hot' if its value of 
exceeds AkeV, 'normal' if > 2keV or 'cool' when 7™J < 2keV. 
Moreover, the subsamples are also separated according to the degree of reg- 
ularity of the gas distribution of its cluster members. A cluster is part of a 
subsample denominated 'quiescent' if the value of n 3 is below the threshold 
value which defines the 25% percentile of the cumulative distribution of the 
power ratios. Similarly, clusters which are members of the 'active' subsam- 
ple have their value of II3 above the threshold which defines the 75% of the 
percentile of the distribution. 
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For a given subsample, the mean spectral temperature profile is defined by 
averaging the profiles of all the clusters of the subsample, with the constraint 
that a profile is part of the mean if there are at least five spectral temper- 
atures with contiguous radial bins which satisfy the constraints of Sect. 3.3. 
This criterion was introduced because it has been found that otherwise there 
would have been clusters with very few annulus bins which would have taken 
part in the construction of the mean of the profile. These clusters with a 
very sparse sampling of their profile lead to a distortion of the average pro- 
file owing to the scale dependency introduced by cooling. If a cluster spectral 
profile is accepted to be part of the mean of a given subsample, then the cor- 
responding cluster emission-weighted temperature profile is also part of the 
average emission-weighted profile for the subsample under consideration. Av- 
erages have been performed by rescaling the cluster temperatures to or 
T e 2 °°[0.5 — 10] = T x , so that it is the rescaled profile which is averaged. This 
allows us to compare consistently mean spectral profiles with mean emission- 
weighted profiles. In order to properly compare with available data profiles, 
the cluster temperature Tx has been calculated by excluding a central region 
of size 50h~ 1 kpc. 

The upper panels of Fig. 6 show the projected spectral temperature profiles as 
a function of r/r2oo- Within a given panel points with different symbols refer 
to mean spectral profiles as obtained from different subsamples previously 
defined; dashed lines are the mean emission-weighted profiles extracted from 
the same subsample. The profiles of quiescent clusters are displayed in the right 
panel, those of the active cluster in the left panel. The top right panel of Fig. 6 
shows that quiescent clusters have scaled profiles which rise toward the cluster 
center, reaching their peak values at r ~ 0.02r2oo and with a steep decline 
thereafter. Spectral profiles follow this behavior, but with a biasing which is 
strongly dependent on the chosen subsample. For hot clusters T$(r) ~ T ew (r) 
whereas for cool clusters this is valid only in the inner regions. Moving outward 
from the peak value the bias can be as high as ~ 20%. The discussions of the 
previous sections suggest that these biasing dependencies can be explained 
as follows: spectroscopic measurements are biased by line emission toward 
lower values than emission-weighted temperatures, this effect is negligible for 
hot clusters but relevant for those clusters of the T < 2 subsample. This 
biasing is strongly suppressed at small radii because of the cooling efficiency 
for these clusters, which has removed most of the cold gas at the cluster 
cores. The same arguments apply to the active clusters of the left panel. Here, 
however, the profiles are shallower than those of the right panel and the peak 
heights are much more modest. Normal and hot clusters have their profiles 
approximately isothermal, in the cluster inner regions spectral temperature 
profiles are now also biased with respect emission-weighted profiles. These 
dependencies of the shape of the profiles on the value of H3 indicate that the 
effects of merging on the gas distribution of the clusters are the main source 
for the differences in the profiles. According to this framework active clusters 
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have profiles much shallower than quiescent clusters because their cores have 
accreted from subclumps a significant amount of cool gas through a number 
of merging events. For small values of the aperture radius (R ap ~ ^2oo/4) the 
differences between the shape of the profiles from active and quiescent clusters 
are not as well defined as those displayed in Fig. 6. This suggests that cools 
gas can significantly accrete into cluster cores only through major merging 
events, where the mass of the subclump is a significant fraction of the cluster 
mass. 

The solid lines in the upper panels of Fig. 6 indicate the 68% confidence limits 
of the best-fit profile of Allen, Schimdt & Fabian (2001a). The profile was orig- 
inally scaled in units of T 250 o and r 2 5oo, here it has been rescaled in units of T 20 o 
and r2oo- F° r the considered sample of clusters < T2500 > / < T 200 >~ 1.31 
and < r 200 > / < ^2500 >— 2.78, with small dispersions. The best-fit profile 
exhibits a decline in the cluster inner region and an approximate isothermal- 
ity between r ~ 0.1r 2 oo and r ~ 0.3r 2 oo- From Fig. 6 it can be seen that the 
simulated profiles are unable to follow this behavior. They have a steep rise 
toward the cluster center, with peak values located at r ~ 0.02r 20 o- Moreover 
they decline outward with radius, whereas the best-fit profile stays approx- 
imately constant. However it must be stressed that the sample of Allen et 
al. is somewhat peculiar in this aspect. Declining temperature profiles have 
been recently measured by a number of authors (Vikhlinin et al. 2005, Pif- 
faretti et al. 2005), using spatially resolved spectra obtained from Chandra or 
XMM-Newton satellites. The rescaled profiles are quite similar, with a decline 
in temperature from its peak value at r ~ 0.1r 20 o toward outer radii. A possi- 
ble explanation for this discrepancy is suggested by analyzing the differences 
between the profiles of different subsamples. From the top-right panel of Fig. 
6 it can be seen that the average profile of cool clusters is steeper than that of 
hot clusters. This dependence of the shape of the profile on the cluster mass 
is a consequence of the scale dependence introduced by cooling. According to 
the results shown in the previous section the rescaled central temperatures are 
expected to be higher for cool clusters than for hot clusters. As a consequence, 
the slopes of the temperature profile at the cluster center will also be higher. 
This dependency of the slopes on the cluster mass is more clearly illustrated in 
Fig. 9 and is observationally confirmed by Fig. 16 of Vikhlinin et al. (2006). In 
this paper the authors analyze the gas and mass-density profiles of a sample 
of 13 low-redshift regular clusters. For this sample spectral temperature pro- 
files were extracted in Vikhlinin et al. (2005), and in the lower panels of Fig. 6 
their best-fit profile is compared against the profiles from the cooling runs (see 
later). Fig. 16 of Vikhlinin et al. (2006) shows that the observed profiles for 
cool clusters are significantly different and steeper than for massive clusters, 
thus confirming the scale dependency found here in the scaled profiles. 

The average profile of hot clusters is the one which is closest to the strip which 
is defined by the best-fit profile of Allen, Schimdt & Fabian (2001a). The agree- 
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Fig. 6. Spectroscopic projected temperature profiles (points) are compared against 
emission- weighted profiles (dashed lines) and best-fit data constraints (solid lines). 
The dashed lines follow from bottom to top the same subsample order of the spec- 
tral points. The profiles are the averages over individual profiles of those clusters 
which are part of a subsample. These have been defined as follows : the notation 
T > 4 means that all the clusters with T^w > 4:keV are part of the subsample, 
while II3 > 75% means that these clusters have a value of II3 above the threshold 
which defines the 75% of the cumulative distribution. The meaning of the other 
notations is similar. The power ratios have been calculated with R ap = r2oo/2. The 
upper panels have profiles extracted from the spectral sample of clusters at redshift 
z = 0.038. The profiles of the lower panels have been averaged over samples at 
redshifts z = 0.116,0.052,0.039 and z = 0.025. The solid lines of the upper panels 
define the 68% c.l. of the rescaled best-fit profile T(r)/T™° of Allen et al. (2001a). 
For the lower panels the solid lines give the 68% uncertainties of the best-fit profile 
to the data of Fig. 18 Vikhlinin et al. (2005); temperatures are in units of Tx, the 
cluster emission- weighted temperature in the energy band [0.5 — lOjfeeV^. 

ment is significantly improved if one considers the profiles of active clusters 
in the left panel. In fact, the average profile of the subsample defined by both 
hot and active clusters is in good agreement with data, with the expection 
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of the two innermost bins. These results suggest that a proper comparison of 
simulated profiles with the best-fit profile of Allen et al. (2001a) must take into 
account the mass and the morphological composition of the data sample. From 
Table 3 (ACDM ) of Allen et al. there are only two clusters with T2500 — QkeV, 
whereas T™° > WkeV for the remaining four clusters. Here the subsample of 
hot and active clusters has 5 members with between 5 and 6keV. This 

implies that the scale dependency of the average profiles because of cooling is 
an important factor to be considered. Therefore in order to perform a careful 
comparison of the simulated profiles with the data of Allen et al. a sample of 
clusters extracted from a simulation volume much larger than the ones used 
here (~ 1200Mpc) is needed. 



Another issue concerns the role of merging events that act on the shaping of 
the cluster temperature profile. The best accord with data is obtained for hot 
clusters which are also active, however the sample of Allen et al. consists of 
six clusters for which lensing mass measurements are in agreement with X- 
ray masses. A result which is indicative of relaxed configurations. It is worth 
noticing, however, that some sort of merging activity must have been at work, 
for at least some clusters of the sample. This is indicated by the shape of 
best-fit profile, which is flat or decreases towards smaller radii. This is at 
variance with the scaled profiles as measured from other samples (Vikhlinin et 
al. 2005, Piffaretti et al. 2005), which have a peak at ~ 0.1r 2 oo and a decline 
outward. A flat profile is however what is found when moving from quiescent 
to active clusters: because of the mixing effects associated with the mergers 
the profiles have their peak height reduced or erased completely in the case of 
hot clusters. This behavior agrees with the findings of De Grandi & Molendi 
(2002), who have investigated the temperature profiles for two samples of 
11 non-cooling flow clusters (NCF) and 10 cooling flow clusters (CF). For 
NCF clusters the authors derive profiles with a core which is approximately 
isothermal, which in the standard scenario is interpreted as a consequence of 
mergers. 



In the lower panels of Fig. 6 the simulated profiles scaled in units of Tx are 
shown. The meaning of the symbols is the same as in the upper panels, the 
profiles have been obtained by averaging over cluster samples at redshifts 
z = 0.116,0.052,0.039 and z = 0.025. The solid lines indicate the scatter 
(~ 20%) in the best-fit profile to the data of Fig. 18 of Vikhlinin et al. (2005), 
which is a sample of 13 low-redshift regular clusters. Spectral fits have been 
performed by Vikhlinin et al. in the energy band [0.6 — lOjfce^, whereas the 
fits here are in the energy band [0.5 — 7}keV. The differences in the spectral fits 
can be considered negligible, since all the simulated clusters have J 1 ™ ^ QkeV . 
The radial dependence of the best-fit profile is given by Eq. 2 of Vikhlinin et 
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al. (2005): 



f 1.07 if 0.035 < r/r x < 0.125 

T/T x = X (18) 

[ 1.22 - 1.2r/r x if 0.125 < r/r x < 0.6 , 

where rx = 1.95/i _1 (Tx/10fcel A ) 1 / 2 . The solid lines have been extrapolated in 
the plots downward of r/r2oo < 0.035 for the sake of a comparison with the 
spectral points. For r/r2oo > 0.1 the profiles are marginally consistent with 
data, at r/r 20 o — 0.2 the simulated profiles are about ~ 10% — 20% lower 
than the central value of the best-fit profile (T/Tx — 0.1). There is a formal 
improvement at r/r2oo < 0.1, where the flatness of the solid lines accounts in 
cool cores for the scatter of the measured profiles. Clearly, a larger sample of 
measured profiles of cooling flow clusters is needed before it would be possible 
to draw statistical meaningful conclusions about the consistency of the profiles 
of simulated cooling clusters with data. From Fig. 15 of Vikhlinin et al. (2005), 
but see also Fig. 2 of Piffaretti et al. (2005), there is a clear concordance that 
the measured scaled profiles reach their peak values at r/r2oo — 0.1, whereas 
the simulated profiles have their peak much closer to the cluster center, at 
approximately r/r 2 oo — 0.02. The disagreement with data can be improved by 
taking advantage of the scale dependencies previously discussed and by con- 
structing a cluster subsample of a given mass and morphological composition. 
However, the general framework predicts a strong correspondence between the 
presence of a cooling flow and a regular cluster morphology (Bauer et al. 2005, 
but for an alternative view see Motl et al. (2004)) 

Because the sample of Vikhlinin et al. consists of 11 relaxed clusters this indi- 
cates that one should compare with data the simulated profiles of the bottom 
right panel, rather than the profiles of active clusters of the left panel, though 
a proper comparison is not possible without a quantitative morphological mea- 
sure of the data clusters. Although appears unlikely that the simulated tem- 
perature profiles of simulated cooling clusters can consistently fit the profiles 
of cooling flow clusters, note that a temperature drop at the cluster center is 
however qualitatively reproduced. 

In Fig. 7 the sums J2i ^Tf are shown as a function of the cluster power ratios, 
each panel is for a different value of the aperture radius R ap . For a given cluster 
the sum J2i ATf is defined as the sum over the radial bins, out to r/r 2 oo = 0.1, 
of AT? = [(T s ( ri ) - T ew {ri))/T™)} 2 . The sample investigated is the T > 2 
in the upper panels of Fig. 6. In each panel the Spearman significance levels 
of the correlations are given. The main result is that there is a significant 
correlation between J2i AXf and II3, a result which indicates that the role of 
merging events in reshaping the temperature profiles can not be considered 
marginal. Moreover, for quiescent clusters the correlation between J2i ATf and 
the cluster dynamical state appears more robust. 
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Fig. 7. The sums ^ AT 2 versus the cluster values II3. For a cluster ^ AT 2 is 
defined as the sum over the radial bins, out to r/r2oo = 0.1, of the squared difference 
AT 2 = [(T s (r) - T ew (r))/T™)] 2 . The clusters are those of the sample T > 2 of 
the upper panels of Fig. 6. Each panel is for different values of R ap with which 
the power ratios are calculated. Open squares (filled triangles) refer to all (25% 
percentile) of the power ratio distribution. The notation Spe= (Spel=) gives the 
Spearman significance levels of the correlations. 



In order to assess the consistency of the simulated profiles with previous re- 
sults, for comparative purposes the work by Loken et al. (2002) was chosen, 
who used an adaptative Eulerian mesh hydrodynamic code. The projected 
emission-weighted profiles at z — 0.038 are shown in Fig. 8; the radial coordi- 
nate is now scaled to rx and the horizontal scale it extends out to r/rx = 0.4. 
The right panel is for the subsample with il 3 below the 25% percentile, left 
panel is for all of the clusters. The solid line indicates the best-fit profile of 
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Loken et al. (2002), 

T/T x = T (l + r/a)- s , 



(19) 



where T = 1.33, a = r vir /1.5 and 5 = 1.6. Their emission-weighted profiles 
were constructed assuming a constant metallicity and using a [1.5 — ll]keV 
bandpass, thus it is appropriate to compare the profile of equation (19) with 
simulated profiles of hot or normal (T > 2) clusters. The profiles of Fig. 8 are 
in good agreement with the best-fit profile (19) over the radial range 0.1 < 
r / r x < 0.4, with r/rx = 0.04 being the lower limit to the fit. The profiles 
of the whole cluster sample are shallower than those of quiescent clusters and 
in better agreement with the profile of Loken et al. (2002). A result which 
indicates that the total sample is more representative of the sample of Loken 
et al. over which the fit was performed. 

A comparison of the simulated spectral temperature profiles presented here 
with the data of Piffaretti et al. (2005) is not directly possible because the au- 
thors have measured their profiles using the cameras on board XMM-Newton. 
Nonetheless an indirect comparison is still possible through their Fig. 6. In 
this figure the authors report, together with their data points, the best-fit of 
equation (19). The results indicate a substantial agreement of the simulated 
profiles with data out to r/rx — 0.4. At small radii (r < 0.05rx) a comparison 
of the simulated profiles with that of Loken et al. (2002) is not possible due to 
the radial range of the fit, though it would have been interesting because of the 
different numerical methods employed in the simulations. If one considers pre- 
vious discussions and the profiles of Fig. 6 and 8 over their entire radial range, 
it turns out that the temperature profiles of cooling clusters is scale dependent 
and can not be considered universal, the profile of hot clusters being shallower 
than that of cool clusters. This is clearly illustrated in Fig. 9 (lower panels), 
where the dimensionless gradients VT = d(T/Tx)/d(r/rx) are evaluated at 
r = 0.05rx for the sample of Fig. 8 and the cluster values are shown against 
Tmw- The correlation of VT against is well defined, the significance levels 
of the linear correlations are found to be above 95%. This scale dependency 
has been found here clearly because of the large size of the numerical sample 
used; this has allowed us to subdivide the sample into subsamples according 
to several cluster properties, while keeping the subsamples with a statistically 
significant number of members. 

Finally, in the top panels of Fig. 9 the dependencies of the mass accretion 
rates M and the central cooling times r c versus the cluster dynamical state 
as indicated by n 3 are investigated. Note that here the power ratios have 
been evaluated at R ap = r 2 oo, while R ap = T200/8 for the power ratios of the 
lower panels. For a given radial bin the mass accretion rate is defined as the 
spherical average of M(r) = Airp g r 2 v r , where v r is the radial velocity of the 
gas. Both M and r c are evaluated at a fixed radius of 50kpc. According to 
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Fig. 8. Average emission-weighted profiles are displayed for different subsamples of 
the simulated clusters at z = 0.038. The notation T > 2 indicates, as in Fig. 6, that 
all the clusters with T™ > 2keV are part of the subsample, while 1I3 < 25% means 
that these clusters have a value of II3 below the threshold which defines the 25% 
of the cumulative distribution. Tx is the emission- weighted cluster temperature in 
the [0.5 - !Q]keV energy band and r x = l.Qhh~ l (T x /IQkeV) 1 / 2 Mpc. Continuous 
line is the best- fit profile of Loken et al. (2002). Left panel is for all of the clusters 
while the right panel is for the cluster subsample H3 < 25%, the power ratios have 
been calculated setting R ap = T2oq/2. 

the standard scenario (Fabian 1994) cool gas at the cluster center will form a 
cooling flow as the cooling time of the gas will be shorter than the Hubble age. 
Spectroscopically the mass deposition rates of cooling flow clusters are found 
to lie in the range ~ 10 — 5OOM yr _1 . In order to mantain a steady state flow 
the inner region of the cluster must remain undisturbed by merging events. 
Therefore one expects to see some degree of correlation existing between the 
central cooling time r c , or the mass deposition rate M, versus a morphological 
measure as n 3 . From the top left panel of Fig. 9 it can be seen that the 
accretion rates M have a widely scattered range of values, without displaying 
any clear dependence on II3. The cause of the scatter is possibly due to the 
adopted definition of M, which requires an estimate of the infall velocity which 
is subject to noise because of the local gas motion. 

The cooling times r c show instead a significant correlation with the values of II3 
and thus with the cluster morphology. The Spearman significance level of the 
correlation is ps = 0.03. For quiescent clusters of the sample the average of the 
cooling times gives < r c > q = 3.75Gyr and for active clusters < r c > a = 43Gyr. 
A Student t-test applied to the two distributions yields p t = 0.08, so that the 
significance level for different means is above 90%. This is an important result 
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Fig. 9. The mass accretion rates and central cooling times of the simulated clusters 
are displayed at z = 0.038 in the upper panels against the power ratios log 10 P^/Pq, 
these have been evaluated at R ap = r2oo- I n the lower panels values at r = 0.05rx 
of the rescaled temperature gradient VT = d{T /Tx) / d{r /rx) of Fig. 8 are plotted 
versus T^®. Left (right) panel is for all (25% percentile) of the power ratio distri- 
bution, with R ap = r2oo/8 for the power ratios. The Spearman significance levels of 
the linear correlation coefficients is below 95% only for the top left panel. 

because it confirms through numerical simulations the association between 
cluster morphologies and the strength of cooling cores. The plot can be com- 
pared with Fig. 8 of Bauer et al. (2005), who have measured from Chandra 
temperature maps the cooling times and the power ratios of 38 X-ray luminous 
clusters. The correlation of cooling times with U 3 is qualitatively reproduced, 
though with several differences. The range of values of the simulated power 
ratios is higher than for data, the latter being £ 10~ 6 . Moreover, the distri- 
bution of the central cooling times is wider for data than for the simulations, 
with r c lOGyr for II3 £ 10~ 7 . The main source of these differences are the 
different aperture radii used by Bauer et al. (R ap = 0.5(/i/0.7) _1 Mpc) and the 
definition of central cooling time which is measured at the cluster innermost 
radial bin, and is distance dependent, whereas here a fixed radius of 50kpc 
was used. 
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A systematic analysis of the evolution of cool cores in galaxy clusters and how 
they are correlated with the cluster merger histories will be investigated in a 
future paper; here the results on the cooling times are used to confirm that the 
temperature profiles of quiescent clusters in Figures 6 and 8 can be identified 
with those predicted by the simulations for cooling flow clusters. 



6 Numerical issues 



In this section it is investigated the dependence of the results presented in 
this paper on the numerical resolution of the simulations. The latter is mainly 
controlled by two numerical input parameters: the number of simulation par- 
ticles and the chosen values for the gravitational softening parameters. The 
most important numerical effect which must be kept under control in numeri- 
cal simulations is the 2-body heating time r r (r), which must be at least of the 
order of the age of the structure and is defined as (Binney & Tremaine 1987) 

TV ~ 6.7 • WGyr ^^——^ ^ md/1QllMQ ^ pd / Pc) i nA > ( 20 ) 



where 0\ is the 1-D dark matter velocity dispersion, G is the gravitational con- 
stant, pd is the dark matter density; In A is the Coulomb logarithm associated 
to the gravitational interaction, and typical values are In A ~ 3; The relax- 
ation time r r can be estimated as follows: for a Navarro, Frenk & White (1995) 
dark matter density profile (Lokas & Mamon 2001) at z = the 1-D central 
velocity dispersion scales as a ~ 950km sec" ^(M^/IO^/i^Mq) 1 / 3 . Moreover, 
with the simulation method described in sect. 2, the dark particle mass rrid of 
the simulated cluster can be expressed as a function of the mass of the cube 
of size L c and thus it scales linearly with M v i r . For a lattice with Nl = 51 3 
grid points it is then easy to find that rrid — 2 ■ 10~ 5 M^ r . Accordingly, eq. 20 
yields: 



The heating time is then only weakly dependent on the cluster mass M vir 
through the dark matter density Pd{r)- For a given cluster the heating time 
decreases as the density gets higher, r r is evaluated at the core radius r c ~ 
0.01r2oo> approximately the resolution limit of our simulations. For the most 
massive cluster of the numerical sample r 2 oo — 2Mpc, e g = 25kpc, In A ~ 2.6 
and pd/ 'p c — 10 5 at the radius r c . The 2-body heating time is then r r ~ 15 
Gyr. Similar values are obtained for the least massive clusters of the sample. 
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Table 3 

Reference values for the four clusters used in the numerical tests. M2oo : cluster mass 
at z = within r2oo in /i _1 Mq, r2oo in units of Mpc and T^® in keV. The power 
ratio H3 is evaluated at R ap = T2oo- e g'- value of the gas softening parameter in 
kpc, N g : number of gas particles. These values are for the high resolution runs. The 
last two columns give the values of the stellar mass fraction f s tar(< r) = M s t ar (< 
r)/M to t(< r) evaluated at r = 0.1r2oo and in units of Qb/£l m , for both the stantard 
(S) and the high high-resolution (H) runs. 



cluster 


M 200 


^200 


7^200 

mui 


n 3 


£ 9 


N 9 


fs(S) 


fs(H) 


003 


6.9 • 10 14 
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0.108 
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6.2 • 10 14 


1.99 


4.6 


-7.32 


17 


212051 


0.104 


0.123 


092 


7 • 10 13 


0.96 


1.29 
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3.3 • 10 13 


0.75 


0.95 
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This shows that, outside of the core radius, the simulations can be considered 
free of 2-body heating. 

Another timescale which is relevant for these simulations is the cooling time 
r c = 3nksT /2A C , where ks is the Boltzmann constant and n is the gas number 
density. According to Steinmetz & White (1997) gas cooling will be affected 
by artificial 2-body heating unless the dark particle mass is smaller than the 
critical value 

M c = 2 • 1O 9 T 6 / O . O5 M , (22) 

where T 6 is the the gas temperature in units of 10 6o J^, /0.05 is the ratio / = 
pg/pd in units of 0.05. For the simulated clusters studied here T 6 lies in the 
range T 6 ~ 10 — 80 and / .os ~ 2 for r £ 0.1r 2 oo- Thus the critical mass M c is 
between ~ 3 • 1O 11 M and ~ 4 • 1O 1O M . This range of values is always above 
the value of m d of the corresponding cluster by a factor ~ 10, so that the gas 
behavior can be considered free from numerical effects. 

In order to support these conclusions, the stability of the final results against 
the numerical resolution of the simulations was tested by running again four 
different cluster simulations, but with a number of particles increased by a 
factor ~ 3. The clusters were chosen in two pairs, with the constraint of being 
the clusters of the first pair among the most massive of the sample and those 
of the second pair among the least massive. The choice of the clusters was 
further dictated by the requirement of being in a relaxed state, therefore only 
those clusters with a value of the power ratio il 3 , evaluated at z = and at the 
aperture radius R ap = r 2 oo, below the threshold value which defines the 50% of 
the cumulative distribution were eligible. The four clusters identified accord- 
ing to these criteria are listed in Table 3, they are identified by their index 
in the sample and have N g ~ 210, 000. The gravitational softening parameter 



37 



LO 

o 



0.1 



003 
008 



0.01 



0.1 



I 

LO 

o 



0.1 



0.01 



144 

092 



0.1 



Fig. 10. Radial dependence at z = of the temperature profiles of the four test 
clusters (see Table 3). Each cluster is labeled by its sample index, the standard 
simulations (sect. 2) are indicated by thin lines and the high resolution runs by 
thick (red) lines. Left panels are for the two most massive clusters and right panels 
for the two least massive. In the two panels are displayed the projected emission 
weighted temperature profiles in the [0.5 — 10] fee V band in units of T^®, the vertical 
lines indicate the value of e g for the simulation. 

for the gas particles is set to 17(10)fcpc for the most (least) massive pair of 
clusters. The other numerical parameters have their values scaled accordingly. 
The stability of the final temperature profiles against the numerical resolution 
of the simulation can be estimated from Fig. 10, in which the projected emis- 
sion weighted temperature profiles in the [0.5 — 10]&eV band are displayed in 
the top panels for the standard resolution and the high resolution runs. For 
all of the test clusters there is a good agreement between the profiles with dif- 
ferent resolutions, indicating that convergence of hydrodynamical variables is 
achieved for N g £ 70, 000. For the less luminous clusters the similarity of the 
temperature profiles between runs with different resolutions shows that the 
peak location of the profiles is adequately resolved in the standard resolution 
runs and is not adversely affected by the value of e g . 

Another numerical issue is the effect of numerical resolution on the final 
amount of cooled gas turned into stars. In the last two columns of Table 
3 for each cluster the values of the star mass fraction f s tar(< r ) = M star (< 
r)/M tot (< r), evaluated at r = 0.1r 20 o, are given in correspondence of the 
standard and high-resolution runs, respectively. The similarity of the values 
of fstar between runs with different resolution suggests that convergence is 
being achieved for f star and that the SN feedback model implemented here 
effectively regulates the amount of cooled gas. These conclusions agree with 
previous findings of V03 (sect. 4.4), in which it was shown the stability of 
fstar against the numerical resolution of the hydrodynamic cluster simulations 
performed there. Finally it is worth noting that cluster X-ray properties are 
largely determined by the high-temperature part of the gas distribution, and 
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are not very sensitive to the amount of cooled gas in the cluster inner regions 
(V03). 



7 Conclusions 

In this paper results from a large set of hydrodynamic/N-body simulations 
of cooling clusters have been used to investigate the relationships between 
spectral temperatures and mass-weighted or emission-weighted temperatures. 
X-ray spectral temperatures have been constructed according to a set of 
procedures in order to reproduce spectroscopic measurements as expected 
from Chandra satellite. Spectral fits have been performed using the single- 
temperature mekal model of XSPEC, with three free parameters: the gas tem- 
perature, the metallicity abundance and the normalization. The main differ- 
ence with respect to previous works which have investigated these issues is in 
the physical modeling of the gas, which is allowed to cool radiatively and is 
subject to feedback heating by SNe. This in turn implies a scale dependency 
of cluster X-ray properties because of the cooling efficiency to convert cold 
gas into stars, which is higher for cool clusters. For instance, the relationships 
between the global spectral temperatures Tg and the mass-weighted temper- 
atures T^ w are found to be significantly biased with respect to those obtained 
from previous runs (ME), in which the modeling of the gas was adiabatic. 
The scale dependency introduced by cooling has a significant impact also on 
the value of T^ w as different cluster overdensities are considered, therefore 
the Tg — T^ w relation is found to be strongly affected also by the value of 
A. Moreover, the scale dependency due to cooling also enters indirectly in 
the determination of the shape of the Tg — T^ w relation because of selection 
effects which limit the sample construction and are dependent on the cluster 
mass. The biasing of spectral temperatures with respect to emission-weighted 
temperatures is found to be similar to that with mass-weighted temperatures, 
but with a wider scatter. 

Another effect which governs the degree of biasing between spectral and mass- 
weighted temperatures is given by merging processes. Spectral temperatures, 
as obtained by fitting single-temperature models, can be significantly biased 
toward lower values than mass-weighted temperatures because of accretion of 
subclumps due to hierarchical clustering which contain cooler gas that mod- 
ifies the photon spectrum. This source of biasing was discussed in detail by 
ME, here it has been analyzed using the power ratios as indicators of the 
cluster substructure. The main advantage of the method is the possibility to 
quantify the degree of substructure on a chosen scale, identified by the aper- 
ture radius R ap . This scale dependency has been investigated by evaluating 
for the considered samples the cluster power ratio n 3 at the aperture radius 
R ap ~ i? A , in correspondence to Tg /T^ w for the cluster under consideration. 
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The results confirm the existence of a strong correlation between the ratio 
Tg /T^ w and the cluster substructure as measured by n 3 . This correlation is 
found with a high significance level for all of the considered overdensities. An 
important result which follows from this analysis is that the biasing associated 
with merging events is independent from the one which follows because of the 
scale dependency due to cooling. Thus, the biasing of spectral temperatures in 
galaxy clusters is governed by two processes, independent of each other, the 
first being due to subclump accretion of cool gas and the second because of 
cooling. 

This two parameter model can explain also the shape of the simulated tem- 
perature profiles. The temperature profile is however not universal even for 
relaxed clusters. This follows because, owing to the scale dependency intro- 
duced by cooling, the central temperature in relative units is higher for cool 
clusters than for hot clusters. This in turn implies for the former steeper 
profiles at the cluster center than for massive clusters. The gradients of the 
rescaled profiles are found to be dependent on T^ w with a high statistical 
significance. This scale dependency of the shape of the profiles on the cluster 
mass is softened or erased entirely if the cluster has undergone a number of 
major merger events. The profiles of these clusters are much shallower than 
those of quiescent clusters. This is clearly due to the effects of remixing of the 
gas due to mergers and these differences have been quantified according to the 
cluster morphology using the power ratios. It has been found that the profiles 
are modified in a significant way only in the presence of major mergers, where 
the mass of the accreting subclump is a large fraction of the cluster mass. 
These results were already suggested by Allen et al. (2001b, cf. sect. 9.2), but 
here they are confirmed through numerical simulations. 

For quiescent clusters the scaled temperature profiles are in good agreement 
with the considered data in the radial range 0.1 £ r/rx ~ 0.4. The profiles 
are not consistent with those measured from a sample of massive clusters ( 
Allen et al. 2001a), but the disagreement can be made less severe if the simu- 
lated sample is constructed by taking into account the scale dependencies due 
to cooling and merging events. A more severe discrepancy is instead found 
at small radii (r/rx ~ 0.1) between the measured profiles and those of the 
simulations. The profiles of the simulations reach their peak values at the 
scaled radial coordinate r/rx — 0.02, which is smaller by a factor ~ 5 of the 
peak location as derived by observations. Moreover, the peak height of the 
simulated temperature profiles is ~ 20% higher than that measured. These 
discrepancies are clearly indicative that the cooling model fails to consistently 
reproduce the observed temperature profiles. In a recent paper Motl et al. 
(2004) argue that hierarchical merging can provide a viable scenario for the 
formation and evolution of cool cores in galaxy clusters. The results of Fig. 
6 indicate that merging events play a key role in reshaping the cluster pro- 
files and support the proposed model. However, it is worth noticing that a 
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significant reshaping of the temperature peak is achieved only in the regime 
of strong merging (II3 £ 75%); moreover this regime has the side effect that 
the central temperature drop is no longer reproduced. 

The measured temperature profiles have a peak location in rescaled radial 
units which is located at a larger distance from the cluster center than in 
the simulations. This is indicative that real clusters, at some stage of their 
evolution, have undergone a significant diminution of their cooling rates in 
their cores. This is in order to provide pressure support and to prevent the 
inflow of gas at a higher entropy toward the cluster center. The suppression 
of the cooling rate is most easily achieved if the gas is heated by some energy 
source. The most plausible candidates are SNe whose feedback energy from 
their explosions can heat the gas to the desired level. However the results 
of the simulations performed here show that this result is not achieved even 
in the case of maximum theoretical efficiency of transfer of energy to the 
ICM. This indicates that the complex nature of the thermal structure of the 
cluster gas is not adequately described by the physical processes incorporated 
in the cooling runs performed here and suggests that additional physics is 
clearly needed in order to prevent overcooling in cluster cores. These results 
agree with previous findings and have prompted many authors to consider 
possible physical sources of additional heating to the ICM. The most popular 
models are (cf. Voit 2005 and references cited therein) energy feedback from 
active galactic nuclei (AGN) and thermal conduction. Because of the existing 
temperature gradients the latter model could in principle provide the necessary 
heating in the cluster central regions (Zakamska & Narayan 2003, Voigt &: 
Fabian 2004), but is currently disfavored by more sophisticated numerical 
simulations which include thermal conduction (Dolag et al. 2004, Pope et 
al. 2005). 

The most favored scenario is thus the AGN model, in which a supermassive 
black hole at the center of the cluster is fueled by the gas of the cooling flow 
and can balance the energy losses of the gas through some energy transfer. 
A popular mechanism involves heating of the gas due to energy dissipation 
of buoyant bubbles, inflated by the AGN near the cluster center and rising 
through the ICM (Churazov et al. 2001, Briiggen & Kaiser 2002, Fabian et 
al. 2002, Ruszkowski & Begelman 2002). The dissipation can occur through a 
number of processes, such as viscous dissipation or turbulence. It is clear that 
incorporating these physical processes into numerical simulations of cluster for- 
mation is a necessary step in order to explain the shape of the observed profiles. 
This is a very demanding computational challenge which has only recently 
begun to be addressed (Briiggen & Kaiser 2002, Briiggen et al. 2002, Bas- 
son & Alexander 2003, Omma et al. 2004, Ruszkowski, Briiggen & Begelman 
2004, Briiggen, Ruszkowski & Hallman 2005, Sijacki & Springel 2006). In 
such a scenario it is difficult to assess how the main findings presented in this 
paper will be modified. Although a proper analysis can only be performed 
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through numerical simulations, several qualitative arguments can nonetheless 
still be used to address the effect of the interaction of the AGN with the ICM 
on the biasing of spectral temperatures. 



In the proposed AGN heating model buoyant bubbles of plasma will rise from 
the cluster center toward the surrounding medium, where the transported en- 
ergy will be transferred into the ICM and converted to heat. This will produce 
some degree of thermalization of the ICM, which in turn implies a reduction in 
the amount of cold gas and thus in the soft X-ray emission, thereby reducing 
the degree of biasing of spectral temperatures . In other words, for a given 
cluster mass scale, any additional source of feedback is expected to reduce the 
biasing. Another issue is if the scale dependency introduced by cooling will be 
preserved in the heating scenario. Some degree of modification will be present, 
since the bubble injection energy, E bu b, is scale dependent itself, according to 
the black hole mass. For example, in one of their models Sijacki & Springel 
(2006) propose E bub oc M^, and therefore in relative terms the energy feed- 
back is higher for massive clusters. The effects on the thermal status of the 
ICM of energy feedback from AGN are however not expected to modify in a 
significant way the scale dependency due to cooling. This argument is mainly 
supported by recent measurements of spectral temperature profiles for a sam- 
ple of cooling flow clusters (Vikhlinin et al. 2005, 2006). The radial behavior 
of the measured profiles for cool clusters is significantly different and steeper 
than for massive clusters, which is in qualitative agreement with what found 
here. These conclusions are valid if the cluster possesses a regular X-ray mor- 
phology, which is indicative of a dynamically relaxed status. In the case of an 
ongoing merger large amount of cold gas could be supplied at the cluster cen- 
ter, thus replenishing the gas reservoir of the black hole and triggering AGN 
activity. In this scenario the thermal status of the ICM surrounding the black 
hole will be due to a complex interplay between the accretion of cold gas and 
the uplifted bubble energy. This suggests that the scale dependency depicted 
in Fig. 4 will be most likely modified in the proximity of a major merger 
event (n 3 — > 0), with a decrease in the spectral bias owing to the feedback 
energy due to accretion and injected into the ICM. Such a complex picture can 
however be addressed self-consistently only through numerical simulations. 



To summarize, it is worth noting, however, that any proposed heating model 
must be able to explain the lack of cold gas at the cluster center as well as the 
observed peak location of the temperature profile. Finally, the results of this 
paper indicate that in order to properly compare X-ray properties of simulated 
clusters with those observed, the mass and the morphological composition of 
the sample must be carefully taken into account. 
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